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We have previously [Phys. Rev. A 65, 043803 (2002)] analyzed adaptive measurements for esti- 
mating the continuously varying phase of a coherent beam, and a broadband squeezed beam. A real 
squeezed beam must have finite photon flux A/" and hence can be signiflcantly squeezed only over 
a limited frequency range. In this paper we analyze adaptive phase measurements of this type for 
a realistic model of a squeezed beam. We show that, provided it is possible to suitably choose the 
parameters of the beam, a mean-square phase uncertainty scaling as {M / k)~^^^ is possible, where k 
is the linewidth of the beam resulting from the fluctuating phase. This is an improvement over the 
(A/'/k)"^/^ scaling found previously for coherent beams. In the experimentally realistic case where 
there is a limit on the maximum squeezing possible, the variance will be reduced below that for 
coherent beams, though the scaling is unchanged. 

PACS numbers: 42.50.Dv,42.50.Lc 
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I. INTRODUCTION 

Optical phase measurements are a valuable means for 
high precision measurement of displacement, for exam- 
ple for gravitational wave detection Two types of 
phase measurement can be distinguished. The first is 
measuring the phase of a beam or pulse of light relative 
to that of a strong local oscillator, which is treated clas- 
sically. The second is measuring the phase shift between 
the light beams in the two arms of an interferometer, 
both of which are treated quantum mechanically. In this 
paper we are concerned solely with the first type. 

Phase measurements are most easily analyzed, and 
thus best understood, for a single mode pulse with a fixed 
phase. The optimal measurement scheme (for all com- 
monly produced states) for such a single-shot measure- 
ment is the so-called canonical measurement 2J. For co- 
herent states, where the quantum noise is independent of 
the quadrature, the canonical phase variance Q asymp- 
totes l/4n in the large n limit, where n is the mean pho- 
ton number. The 1/h scaling is referred to as the stan- 
dard quantum limit (SQL). On the other hand, for more 
general states where the noise is quadrature-dependent, 
the canonical phase variance may approach the Heisen- 
berg limit 1.89/^2 3]. 

Unfortunately it is not possible to achieve canonical 
measurements with linear optical elements (unless one 
discards most of the measurement results Q). One so- 
lution that has been developed is to use feedback. The 
field is combined with a local oscillator at a beam splitter, 
and partial results during the measurement are fed back 
to adjust the phase of the local oscillator for measuring 
the field in the next part of the pulse. This approach has 
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been extensively studied for single-shot measurements, 
and it has been shown that it can achieve scalings al- 
most at the Heisenberg limit &]. It has also been ex- 
perimentally realized Q, verifying an improvement over 
non-adaptive (heterodyne) detection. However, because 
the experiment used coherent states, it was not possible 
to verify a scaling better than the SQL. 

In practice it is easier to produce a continuous squeezed 
beam rather than a squeezed pulse of light. This mo- 
tivates considering continuous, rather than single-shot, 
measurements. In the continuous case, if the phase to 
be measured is constant, then the variance will become 
arbitrarily small with time. To obtain a nontrivial result, 
it is necessary to consider a phase which varies in time. 
We wish to determine how accurately the measurement 
scheme estimates this varying phase. The simplest model 
for a varying phase is a Wiener process; that is, its rate 
of change is Gaussian white noise with intensity k, giving 
rise to a Lorentzian lineshape of the beam with width k. 

Continuous measurements of this type were considered 
m Refs. 0,13. In Ref. 0, coherent beams and broadband 
squeezed beams were considered, and a simple method for 
filtering the data to obtain the phase estimate was used. 
A scaliirg law for the optimal variance (or mean-square 
error) of the squeezed beam was derived which showed an 
improvement over the coherent beam result. However, as 
was pointed out in Ref. jSj], this analysis had two short- 
comings (which tend to counteract each-other). First, it 
considered only the photon flux due to the coherent com- 
ponent of the beam, whereas in fact the photon flux from 
the broadband squeezing is strictly infinite. Second, its 
filtering technique ignored the fact that phase informa- 
tion may be obtained from the photocurrent noise j^. 
For broadband squeezing the information from the noise 
is strictly infinite, and would allow the phase to be de- 
termined exactly, modulo n. 

Ref. adopted a more sophisticated (Bayesian) filter- 
ing of the state, but restricted its attention to coherent 
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states. It found different results for Af ^ k, though the 
asymptotic variance of J_1/2)a/k/A/' for A/" ^ k was the 
same as found in Ref. [3]. An asymptotic improvement 
by a factor of 1/V2 over non-adaptive (heterodyne) de- 
tection was also confirmed. 

In this paper we apply a Bayesian approach to phase 
measurements on a narrowband squeezed beam. This es- 
timation procedure optimally uses the phase information 
obtained from both the mean field and from the noise. 
The narrowband squeezing ensures that the photon flux, 
and the rate of information accumulation in the pho- 
tocurrent record, are finite. To obtain analytical results, 
we concentrate on the asymptotic regime. We find an ap- 
proximate analytical expression for the optimal scaling of 
the phase variance as (A/'/k)^/^, and confirm this via nu- 
merical simulations. Our analysis confirms that there is 
a scaling advantage over the SQL, although it is slightly 
less than we had previously thought [31 ■ 

The remainder of this paper is organized as follows. In 
Sec. m we give a simple explanation of the scaling, and 
explain why it is different from that in Ref. 7] . Then we 
give the Bayesian analysis of measurements with feed- 
back in Sec, mil and of heterodyne measurements in Sec. 
IIVI We revisit the scaling based on the Bayesian treat- 
ment in Sec. then give numerical results in Sec. IVII 
We conclude in Sec. IVIII and give additional details of 
the derivations in the appendices. 



II. SIMPLE EXPLANATION OF SCALING 

In this section we give a simple explanation to predict 
the scaling of the phase variance under continuous adap- 
tive measurements. First we discuss the simple case of a 
coherent beam, then we proceed to the case of squeezing. 



A. Coherent states 

The configuration considered is as in Fig. ^ although 
for the moment everything to the left of the dashed line 
should be ignored. The coherent signal beam has photon 
flux (i.e. mean number of photons per unit time) M. A 
phase shift 9{t) (to be estimated by the experimenter) 
is imposed on the signal, and a known phase shift of 
$(t) may be imposed by the experimenter on the local 
oscillator. These are combined at a beam splitter, and 
the difference photocurrent I{t) gives a measurement of 
a quadrature of the signal. Feedback may be used to 
determine the phase <I>(t) for the local oscillator using a 
signal processor to filter the signal I{t). The estimate of 
d{t) is also determined based on I{t). 

For fixed system phase and a coherent beam, the vari- 
ance for adaptive measurements over time interval At is 
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FIG. 1: The model of the experiment. The phase shift 9 is 
imposed on the continuous beam, and a phase shift <E> is im- 
posed on the local oscillator. These beams are incident on 
a 50:50 beam splitter, and the difference photocurrent I{t) is 
determined. The processor then adjusts $ based on I{t). The 
simplest case is when the signal beam is in a coherent state, 
but for more accurate phase estimation a squeezed beam, pro- 
duced by the apparatus to the left of the dashed vertical line, 
is used. This consists of a parametric down-converter char- 
acterized by intensity damping rate 7 and x''^' nonlinearity 
parameterized by r. 



This is simply the standard result for adaptive measure- 
ments on a coherent field with a mean photon number 
of AW 13 . If there is an estimate of the phase at time 
t, then the phase estimate at time t -\- At may be taken 
to be a weighted average of the estimate at time t and 
the estimate obtained from data in the time interval At. 
The variance at time t + At will then satisfy 
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l/af + 4AtAf' 



(2.2) 



Now we assume that the phase fluctuations on the sig- 
nal are white: 



(2.3) 



where C is Gaussian white noise satisfying {C{t)({t')) — 
6{t — t'). Taking account of these fluctuations, we expect 
the variance at time t + At to be 
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't+At 



l/fj^ + AAtN 



nAt. 



(2.4) 



Provided nAt < al < l/(4At7V), the total change in 
the inverse variance is 



The steady-state variance is therefore 



(2.5) 



(2.6) 



AAtN' 



(2.1) 



This is just what was obtained in Ref. 0- 

In the case of heterodyne measurements, the variance 
for measurements over time interval At is twice that in 
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Eq. 1111) Then the last term in Eq. (ESJ is 2AtAf, 
and the steady state variance is 



V2 



(2.7) 



Thus the adaptive measurements give a l/-\/2 reduction 
in the variance over heterodyne measurements. 



B. Squeezed states 

We may use a similar method in the case of squeezed 
states. In this case we have a number of new features to 
the model, indicated in Fig. 1 by the apparatus to the left 
of the dashed line. The squeezed beam is produced by a 
cavity with decay constant 7. The beam has squeezing 
parameter r and coherent amplitude E. The flux is given 
by [13 



— sinh^ r. 



(2.8) 



For squeezing which is not too large, the phase variance 
for adaptive measurements over time interval At is, using 
the simple estimation technique of Ref . , 

(2.9) 
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That is, the variance is reduced by a factor of e^^ from 
what would be obtained for a coherent state. The equiv- 
alent of Eq. I|2.5|l is then 



-At—r + e^^'AtE' 



(2.10) 



Provided the dominant contribution to the flux is the 
coherent component, the steady-state variance is 
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(2.11) 



There is a limit on how large the squeezing can be 
before this approximation fails. In Ref. 7J the limi- 
tation considered was the error in the feedback phase, 
which causes the measurement to not be exactly on the 
squeezed quadrature. This limitation gives an overall 
scaling of {k/NY^^. 

Here we consider the additional limitation due to the fi- 
nite squeezing bandwidth. For measurements performed 
over a time scale shorter than I/7, no squeezing will be 
observed. In fact, it is necessary to perform measure- 
ments over a time scale of order before squeezing is 
observed. For large, this time is approximately the 
reciprocal of the decay constant for one of the cavity 
quadratures [see Eq. H3.2|) in the following section] . How- 
ever, it is not obvious why this decay constant is impor- 
tant, because it relates to the antisqueezed quadrature. 
We give a more thorough explanation of this time scale 
in the following section. 



The phase information used to obtain the phase esti- 
mate will be from a finite time interval. For times before 
t — ct^/k, the system phase will differ from the current 
system phase by an amount comparable with the phase 
uncertainty. Therefore, the majority of the phase infor- 
mation must be taken from a time interval of length ap- 
proximately ct^/k. We use x to denote the inverse of this 
time interval, so x /c^ . 

In order to observe squeezing, this time interval must 
be longer than the time scale e''/7, so we require ct^/k > 
e'"/7. Here we are concerned with scaling, so we ignore 
multiplicative constants in the rest of this section. In 
order to obtain the scaling we also need to consider the 
limitation due to the squeezing contribution to the flux. 
From Eq. H2.8|l . we have M > 76^'', so 
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(2.12) 



In order to obtain the minimum phase variance, we need 
to take the squeezing r to be as large as possible con- 
sistent with the Eqs. (|2.11l) and 12.12|l . That gives 
~ (AA/k)1/8, so 



K \5/8 
Jf) ■ 



(2.13) 



Note that if we had used the time scale I/7, rather than 
e''/7, then we would have obtained the inequality ct^/k > 
1/7. Then Eq. H2.12|l would become > e'^^'n/Af, and 
the maximum squeezing would be e'' - (TV/k)!/^. rp^^jg 

then would give the variance scaling as ~ (k/A/')^^'^, 
which is identical to that in Ref. 0. Thus we can see 
that the extra is essential to obtaining the different 
scaling here. This scaling will be shown to be correct 
numerically in Sec. IVII 

In practice it is not possible to achieve arbitrary 
squeezing; typically the maximum value of e^'' achieved 
is about 2 [l4|. In that case, we simply have the variance 
given in Eq. 1)2. 11|) , with r taken to be the maximum ex- 
perimentally achievable value. That is, the scaling is the 
same as for coherent states, but the variance is reduced 
by a factor of e"''. 

Note that we have the limitations on 7 (ignoring con- 
stant factors) 



„2r 



-<1< e---. 



K 



(2.14) 



These limitations may be satisfied provided Af/n > e^^; 
for e^*" ~ 2 this limit is about 16. Experimentally there 
is only limited control over the value of 7. However, for 
larger Af /k there is a wide range of values for which a 
reduction in the phase variance should be observed. 



III. ADAPTIVE MEASUREMENTS 

Now we give the detailed Bayesian analysis of the phase 
estimates for narrowband squeezing. We will use this to 
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justify some of the steps used in the previous section, as 
well as to estimate the additional phase information that 
may be obtained from the noise. A continuous squeezed 
beam produced by a cavity may be modeled by the op- 
erator equations in the Heisenberg picture 



^ - -i7(l + £)/2 + 
^ = -y7(l - £)/2 + 



(3.1) 
(3.2) 



/ = cos(<I> - 6*) [^x - I] + sin($ - 0) [^y + E - ff\. 

(3.3) 

These equations are equivalent to Eqs. (4.48) and (4.54) 
in Ref. The quantities x and y are the two quadra- 
tures of the cavity field, ^ and fj are the quadrature noise 
operators, 7 is the cavity decay constant (equivalent to ki 
in J ^^'^ e is a constant related to the usual squeezing 
parameter by 
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1 -£ 



(3.4) 



/ is the output quadrature at angle <^ — 6] the photocur- 
rent measured corresponds to the measured value of this 
operator. Note that the squeezed quadrature here is x 
(rather than y, as in 12]). Also we have added the 
displacement E, so we obtain a squeezed coherent field 
rather than a squeezed vacuum. 

These equations may be solved by using the Wigner 
distribution, and replacing the quadrature operators with 
the corresponding quadrature variables for the Wigner 
distribution. Also the output quadrature / is replaced 
with the detected photocurrent /. 



i = -a;7(H- £)/2 + 

ij = -y7(i - £)/2 + y/iv 

/ = cos($-0)[V7a;-^]- 



sin(<i> 



(3.5) 

(3.6) 

e)[^y + E-ii. 

(3.7) 



Here ^ and 77 are Gaussian increments satisfying 

mm)^mri{t'))^5{t-t'). 

In order to apply the Bayesian approach, we assume 
for the moment that both $ and 9 are constant. The 
results obtained using these assumptions should also be 
accurate for cases where these phases only change by a 
small amount over the time interval considered. We dis- 
cretize the equations to give 



Ax = -xAt7(l + e)/2 + ^\hfAt 



Ay = -yAt-/{l - e)/2 + fj^/^/At 
I At = {cos($ - 6')[V7Atx - I] 



(3.8) 
(3.9) 



+ sm{(S> ~ e)[^/^{y + E') ~ fi]}^/Ai. (3.10) 

Here ^E' = E, and we have used tildes to indicate 
that the stochastic increments have been replaced with 
Gaussian random variables (with mean and variance 

!)■ 



At all times the experimenter's knowledge can be rep- 
resented by a probability distribution for x, y, and 9 
which is Gaussian for x and y: 



P{x,9) = P{9) 



^dct G{9) 1 [g_g(g)]Tg(g-)jg_g|-g)| 



27r 

P{9)P{x\9). 



(3.11) 



where x = {x,y)^ and x{9) is the mean of x. P{9) gives 
the correct probability for 9 averaging over x and y, and 
the inverse of G{9) is the covariance matrix for x for a 
given 9. 



A. Phase information 

First we consider the update to the probability dis- 
tribution due to the information from the measurement. 
At each time step we update the probability distribution 
using Bayes' rule 8]: 



P{x,9) P{x,9) X P{I\x,9). 



(3.12) 



Here the constant factor 1/P{I) is omitted, because we 
can normalize at the end of the calculation. The proba- 
bility P{I\x, 9) is given by 



P{I\x,9) 



■ exp[-Ai{/ - ^[cos($ - 6')a 
-sin{<^-9){y + E')]Y/2\ 



■ exp 



^l^{x~BfAA^{x^B) 



where 



A = 



B = 



AI 



1/c 
l/s 



(3.13) 



(3.14) 



where c = cos($ — 0), 
We therefore obtain 



2V7 

s = sin($ - 9) and A/ = / - sE. 



P{I\x, 9)P{x, 9) cx VdetG 

Xexp|-i[(f-i(2))Tg(2)(^_ -(2)) 



-(5(2))Tg(2)5;(2) + ^Tq- ^ ^i^tB^AA^B]] , (3.15) 



where 



G(2) = G + -lAtAA^ 



(G(2))-\Gx-K7AiAA^B), 



(3.16) 
(3.17) 



Hence the updated probability distribution for 9 is 

1, 



detG 



■ exp 



^-{x^Gx 



detG(2) 

+7At(S^A)2 - (x(2))^G(2)x(2)]| . (3.18) 
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In the limit of small At we obtain the differential equa- 
tions 

dG = -fdtAA'^, (3.19) 
dx = dtG'^[y^AIA~-/AA^x], (3.20) 



dt , 



#ogP(0)] = if - y (A/ ~ V7^^5;)2 



(3.21) 



where if is a constant; K may be ignored, as it only 
changes the normalization. 



B. Increments in x and y 

Now we take account of the increments in x and y. 
Given the measurement result /, we have the restriction 
that 



c^ + sfi= VAt[y/j{cx + sy) - A/]. 
We define the new variable 

/i = -(^ - cVM[y/jicx + sy) - AI])/s 
= {fl- s\/~Ki[^{cx + sy) - A/])/c. 



(3.22) 



(3.23) 



This is simply a Gaussian random variable with mean 
zero and variance 1. In terms of this the new difference 
equations are 

Ax = -xAt-i{l + £)/2 - cAi[^A/ - -f{cx + sy)] 

- sfx^fiKt (3.24) 
Ay = -yAt^{l - £)/2 - sAt[^A/ - ^{cx + sy)] 

+ cfi^f^t. (3.25) 

The deterministic part of the increment gives the map- 
ping X ^ Cx — ^AtAIA with 

C=t- At7[(l + £S)/2 - AA^l 
1 



-1 



(3.26) 



The stochastic part increases the covariance matrix for x 
according to G i-^ (G^^ + 7Affl)"^, where 



n 



-sc 



-sc 

^2 



(3.27) 



Overall, we update the covariance matrix and mean to 
G(3) = [G(G(2))-iG + 7Atn]-\ 



?(3) 



Cx 



(2) 



^AtAIA. 



(3.28) 



In the limit of infinitesimal At we again obtain differen- 
tials, which when added to those at the end of Sec. Ill A 
give 

dG = -idt[AA^ + G + £(I]G + GS) /2 

^{AA^G + GAA^ + GI{G)l (3.29) 
dx = -■^dt{t + eT.)x/2 

+ ^dt{G'^ -t)A{AI - ^A^x). (3.30) 



C. Solution 

From Eq. (|3.21() . the probability distribution for the 
phase is obtained by integrating over (A/ — ^A^x)"^. 
From the reasoning given in Appendix ^ we may es- 
timate the final phase variance by determining the ex- 
pectation value for a given phase Q, and expanding to 
second order in 9 — Q. The probability distribution for / 
is determined based on Q, and we determine the proba- 
bility distribution for 9 based on the measurement result 
/. We interpret O as the actual system phase, and 9 as 
a dummy variable used for the probability distribution. 

For fixed system phase and feedback phase, G will 
reach an equilibrium value. To estimate this equilibrium 
value, we rotate x via 



R = 



c s 
—s c 



(3.31) 



which is defined to simplify AA^ and 11. Using bars to 
denote the rotated values of variables, so for example 
G"^ = RG~^R\ we have from Eq. H3.29|l at steady state 



G-i 



1 




G-^ + G-^ + £(EG^i + G-il])/2 



" 1 


■ 


G-1- 


l-G-i 


■ 1 


■ 


+ 


■ 


" 



















1 



with 



~2sc 



(3.32) 



(3.33) 



Using the notation a = (G ^)oo, b — {G ^)oi, d — 
(G"i)ii, X = {1 - e)/2 + s^e, Y = see, we have the 
simultaneous equations 



= 2{aX + bY), 
b^ + 2{Xd-Yb) = 1, 
ba = Y{a + d). 



Solving for gives 



= 2[X^ + Y^ + v/(A2 + r2)2+y2] 



(3.34) 
(3.35) 
(3.36) 



(3.37) 



Using that equilibrium value, the solution for ^j^A^x 
is (for large time, so initial conditions may be ignored) 

^A^x^-lA^ I exp{7[(l + £S)/2 
Jo 

+ {G-^ - l)AA^]{u-t)) (G^^ - t)AAI{u)du. 
= 7Re [vi due^^("-*)A/(M)^ , (3.38) 

where 

r2 = a-l + i(a-2A:)/A, (3.39) 
2A = a-fiA, (3.40) 
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with AV2 = v/(^^TT2)2Tr2 _ + r2). 

Expanding (AI—y/jA^x)'^ and taking the expectation 
value gives 

((A/-V7^^i)2) = (A/2(i)) 

- 27Re (^nj^ dwe'^^("-*)(A/(i)A/(?i))^ 



^0 



(3.41) 

Determining {I(u)I{v)) gives 

= (1 - 52){7(a;(^.)x(z;)) - VlMu)^{v)) 

+ (y{v)7j{u))] + {fj{uUv))}, (3.42) 
where S = sin((f> — 0). The solutions for x and y are 

x{t) = 77/* e^(^+^)("-*)/2^dv, (3.43) 



y(t) = / e''^^-'^^''-''>^^7]dv. (3.44) 
"'0 



Using this we obtain 



2^2 _^ c^2 jy^p-7(l-e)|«-t,|/2 
1 

+ (^2 _ ^^_2£_g-7(i+c)l«-«l/2^ (3 



{Iiu)I{v)) = (5(u - w) + S^E^ + S'-— — e 



and 



(A/(u)A/(v)) = S{u ~v) + iS- sfE'^ 
1 - e 

(3.46) 

Using this result gives 

((A/ - ^A^xf) = 5(0) + 7/(r, $ - 6, - 6) 



[l-Rc{n/A)]\S^sfE^, 



(3.47) 



where /(r, $ — 9, — 9) is a function independent of 
7 and E; the explicit form will be given below. The 
term on the second line of Eq. H3.47|l may be identified 
with the phase information obtained from the mean field 
(as it is proportional to i?^), and is the most important 
contribution. Expanding to second order, {S — s)^ ~ 
cos2($-9)(6i-9)2, while 



[1 - Re(r!/A)]^ = (c'e-^^ + s^e^^)-' 



(3.48) 
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FIG. 2: The maximum of (?(r, <1> — 9) divided by e^^ as 



function of r. 



In this expression we can replace s with 5, because we 
only need this term to zeroth order in 9 — Q. We then 
have 

((A/ - ViA^xf) ^ „ e)^ 

+ 7/(r, $-9,0-9) + (5(0) + 0{{9 - 9)^), (3.49) 

where S = (1 - S'^)e~^'- + S'^e'^\ 

The additional term 7/(r, $ — 9, — 9) gives the phase 
information obtained from the noise. Its explicit form is 

/(.,a>-9,.-9)^Re LL + _ 



1-e 



+{S^ - 1) 



1 



1 



A + (l-e)/2 
r2V2A+ ir^p/a - 2VL 



(3.50) 



A+(l + e)/2 

Expanding this function to second order in 9 — Q gives 
/(r, $ - 9, e - 9) = const. + .g(r, $ - 9)(6' - 9)^ 



O((0-9)3), 



(3.51) 



where 



2sinh^(r/2) , . , 

.g(r, $ - 9) = ^ ^ [(1 + cosh r) A| 



+2sin-'($ - 9) cosh r 1 



coshr 



(3.52) 



and ae and Ae are equal to a and A for 9 — Q. Numerics 
indicate that g(r, $ - 9) < e^''/4 (see Fig. therefore 
this term may be ignored provided ^e^^ is small com- 
pared to E'^ I'R. 

In Appendix^ the average over 9 is also taken. We 
take $ — 9 to be a constant, rather than taking $ to be 
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independent of 8. To obtain \ogP{9), we take the time 
integral of (A/ — y/jA^x)"^ times —1/2. Therefore, for 
measurement over a time interval At we have 



cos2($ - 6) 



Ata^ 



+ jg{r,^-Q). 



(3.53) 



The first term gives the phase information due to the 
coherent amplitude, whereas jgir) gives the phase infor- 
mation due to the noise. Note that lim,,^o5(^) 0, so 
the phase information from the noise is zero for coherent 
states, as we expect. 

An essential point is the time required for the system to 
come to equilibrium. The real part of the time constant 
7A is a7/2. For small s, we have a/2 « X, so Re(A) w 
l/{e^ + 1). Therefore, the approximations made will only 
be accurate provided the measurement is made over a 
time period that is long compared to {e^ + Ignoring 
the constant 1, this is the e'^/j time scale used in the 
previous section. 



IV. HETERODYNE MEASUREMENTS 

We may derive similar results for the heterodyne case. 
The corresponding equations are 



X = ~xj{l + e)/2 - 
y = -yj{l - e)/2 + /^(??*'^ + V' 
h = cos($ - e)[y^x - ^(1)] 

+ sin($ - e)[v/772(y + E') - 77(1)] 
/2 = - sin($ - e)[y^x - e^^^] 

+ cos($ - e)[^fi/2{y + E') - 77(2)]. 



(4.1) 
(4.2) 

(4.3) 

(4.4) 



The quantities Ii and I2 are the two Fourier components 
of the photocurrents at the frequency at which the local 
osciUator is detuned from the system Because of 

that detuning, $ here is arbitrary. 

To see that the results do not depend on the value of 
$, note that we may represent the measurement results 
by the complex current / — e**(/i + 1/2) given by 



e*«{v/7/2[a; + i{y + E')] -v^- v^}. 



(4.5) 



where v\ = 
g2*(e-*)^^(i) 



(^(1) + j^(2) + ^(2) + i^(i))/2 and V2 



(2). 



irj 

-^(2)_ 



iri^^^)/2 are independent com- 



plex Gaussian random variables. Note also that 
i + = -X7(l + e)/2 - iy-f{l - e)/2 + V^i'i- (4.6) 

Following the derivation in the same way as above 
yields (see Appendix IB|| 



E^ 



Ata^ 1 + e- 



-2r 



27/i(r), 



(4.7) 



As in the adaptive case, the first term gives the phase 
information due to the coherent amplitude, and 2jh(r) 
gives the phase information due to the noise. Again, 
limr^o h{r) — 0, so the phase information from the noise 
is zero for coherent states. For large r, h{r) scales as e^; 
thus the phase information from the noise will be negligi- 
ble unless the contribution to the photon flux due to the 
squeezing is dominant. In the derivation in Appendix IbI 
the (e** + l)/7 time constant again appears, indicating 
that this is the appropriate time scale to observe squeez- 
ing in this case also. 



V. SCALING REVISITED 

We now use these results to give a more detailed deriva- 
tion of the scaling for the phase variance. The derivations 
given in Sees. Illll and llVI are based on fixed system phase 
and feedback phase (for the adaptive case). However, 
provided these phases vary by a small amount over some 
time interval At > (e''-f l)/7, these equations should still 
give an accurate estimate of the phase information over 
this time interval. 

First we consider the adaptive case. Provided the feed- 
back phase is sufficiently accurate, the cos in Eq. (|3.53() 
may be replaced with 1, and S « e"^'". In addition, pro- 
vided ^e^ <C E^ (which will be true asymptotically unless 
the majority of the photon flux is from the squeezing), the 
second term which gives the phase information from the 
noise may be ignored. For continuous measurements, Eq. 
(|3.53|) gives the change in the inverse variance. Therefore 
we have 



1 



1 



(5.1) 



Taking into account the system phase varying as = 
^/k(, we obtain the additional term — k/ct* given in Eq. 

arm . 

We may then follow the derivation given in Sec. HTIto 
obtain the scaling H2.13|l . There are two assumptions 
used in this derivation which we should justify. One as- 
sumption is that the phase information obtained from 
the noise is negligible. For the scaling obtained in Sec. 
ni £;VS - A/'e^'-, whereas je^'' < Me''. Thus the as- 
sumption that little phase information is obtained from 
the noise is reasonable. 

The other assumption is that S ~ e"^''. To check this 
condition, note that sin2($ — 6) may be made as small 
as order cr^ via feedback. Therefore 



sin2($ - 9)6^ 



and 



1/4 



3/8 



where 



(-^) >sin2(<i>-e)e2 



(5.2) 



(5.3) 



h{r) =coshr- (e2 + l)-i/2. 



(4.8) Thus ^ e , as required 
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For the case of heterodyne measurements, the equiva- 
lent of Eq. ifTTUIl is 



1 

At' 

Solving for gives 



1 



-2r ■ 



-2r 



(5.4) 



(5.5) 



Here we have used the approximation that JV « £^^/4, 
which simply means that there is little photon flux from 
the squeezing. If there is no limit to the squeezing, the 
asymptotic variance should be 



VI. NUMERICAL RESULTS 



(5.6) 



As it is predicted that the phase information from the 
noise should be negligible it is reasonable to consider a 
non-Bayesian method which does not use phase informa- 
tion from the noise. The phase estimation method used 
for the majority of the calculations was similar to that in 
Ref. 0- The quantities At and Bt were calculated as 

At= f eX("-*)e**/(u)dM, (6.1) 

^ — OO 

(6.2) 



Bt = - 



gX(i'-t)g2i*_ 



Note that these quantities are not related to the A and 
B used in Sec. 11111 The x used here is equivalent to 
that used in Sec. IllBI in that the time scale over which 
previous measurement results are used is l/x- 

A good phase estimate may be obtained as argCt, 
where Ct = At + X^tA^ 0- It was found that poor 
results were obtained if arg Ct was used in the feedback, 
so the feedback used was 



m = arg(Ci-^A*^ 



(6.3) 



for a suitably chosen 6. 

In the calculations it is not necessary to independently 
vary Af and k. We may scale the time by n, so we obtain 
the dimensionless parameters Af/n, j/k and x/k. The 
other parameters which we may vary are r and S; these 
are already dimensionless. We may predict scalings in 
terms of these dimensionless parameters: 



cr 

7 

K 



e 
X 

K 



(6.4) 



The scaling of x/k is predicted from the fact that we use 
phase information from a time interval ~ ct^/k. 
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FIG. 3: The optimal values of various quantities as well as the 
variance for adaptive measurements with arbitrary squeez- 
ing. The values of a'^ / {M / k)^^^ are shown as the solid line, 
e'' /{H/tiY'^ is shown as the dotted line, {-) / k) / {H / K.f^'^ is 
shown as the dashed line, (x/«;)/(A/'/fi;)^''* is shown as the 
dash-dotted line and 5 / (N / nY^^ is shown as the crosses. 



Numerical calculations were performed for a range of 
values of M j n,. For each value six alternative measure- 
ment schemes were considered: 

1. Adaptive measurements with arbitrary squeezing; 

2. Adaptive measurements with limited squeezing; 

3. Adaptive measurements on coherent states; 

4. Heterodyne measurements with arbitrary squeezing; 

5. Heterodyne measurements with limited squeezing; 

6. Heterodyne measurements on coherent states. 

In each case the appropriate parameters were adjusted 
to minimize the variance. For all cases it was necessary 
to optimize over xl for the adaptive measurement it 
was also necessary to optimize over b. For the cases with 
squeezing the values of r and 7/K were also optimized 
over. In the cases with limited squeezing there was the 
additional restriction that e^"^ < 2. This limit corre- 
sponds to the typical maximum squeezing which may be 
achieved in the laboratory. 

Results were obtained by integrating the system over 
a time interval of 130/x, and determining the average 
phase variance for the data from time 30/x to 130/x. 
A total of 2^*^ independent integrations were performed, 
so the total number of effectively independent samples, 
including those from the different integrations as well as 
those from different times within the integration, was ap- 
proximately 10^. 

The results for adaptive measurements on arbitrarily 
squeezed states are shown in Fig. O In this figure the 
various quantities were divided by the predicted scalings, 
in order to check that these scalings are correct. In each 
case we find that this ratio is of order 1, so the predicted 
scalings are correct. We also obtain the scaling for S as 



{Af/K 



,U/4 



(6.5) 
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n/k 

FIG. 4: The ratio of the value of K/a^ to x (continuous line), 
and the ratio of 'y/e^ to x (dashed line). 



It does not appear to be possible to predict this scaling 
in the same way as the scalings for the other parameters 
given here. 

The above theory predicts that x should be approxi- 
mately equal to k/(T^ in order to minimize the variance. 
Also, the time scale 1/x should be larger than the 
time scale required to obtain squeezing. To check these 
predictions, the ratio of k/ct^ to Xj ^ 'v^^ll as the ratio 
of 7/e'' to X are shown in Fig. 21 From these results, 
the optimal value of x is near k/ct'^, but there can be as 
much as a factor of two difference. We find that 7/e'' is 
larger than x^ ^ predicted. However, the difference is 
not great, and in some cases 7/e'' is only slightly larger 
than X- 

The phase variances for each of the six measurement 
schemes are shown in Fig. |S| In this figure the variances 
have been multiplied by M / k in order to more clearly 
show the scaling constants. In each case except for adap- 
tive measurements with arbitrarily squeezed states, it is 
clear that the variance is scaling as n/N . 

It is str aightforward to predict the asymptotic values of 
(7^ N / K, under each of the measurement schemes using 
the results given in the preceding sections. These predic- 
tions are given in Table HI together with the asymptotic 
values estimated based on the numerical results. In each 
case, the predicted and calculated asymptotic values are 
within 1%. 

Another issue, raised at the end of Sec. II B, is the 
range of possible 7 values for which the benefit from 
squeezing may be observed. In that section we derived 
the upper and lower bounds on 7 in Eq. (|2.14(l . The 
range of values of 7/K such that the phase variance was 
still within 10% of its minimum value was calculated nu- 
merically for adaptive phase measurements, and is shown 
in Fig.|Sl (For the first two data points there was no lower 
bound found numerically, because the variance was close 
to that for a coherent state.) 



1.6f 
1.4 



1.2 




FIG. 5: The phase variance multiplied by ^Af/n. The results 
for adaptive measurements on arbitrarily squeezed states, 
states limited to e^"^ < 2, and coherent states are shown as 
the solid, dotted and dashed lines, respectively. The results 
for heterodyne measurements on arbitrarily squeezed states, 
states limited to e^"" < 2, and coherent states are shown as 
the crosses, plusses and asterisks, respectively. 




FIG. 6: The range of values of 7/K for which it is possible 
to obtain phase variances within 10% of the minimum val- 
ues with limited squeezing (e^"" < 2) and adaptive measure- 
ments. The ranges are shown by the vertical solid lines be- 
tween the crosses, and the predicted upper and lower bounds 
2A/'/(Ksinh^ r) and e^'~ ^/JTJk from Eq. H2.14|l are shown by 
the diagonal solid lines. 

The scaling of the upper and lower bounds found nu- 
merically is the same as the analytic bounds, though the 
range is slightly smaller. For the largest value oi N / n 
tested, there is a range of nine orders of magnitude for 
7. This demonstrates that, even if it is not possible to 
adjust 7 experimentally, there will be a wide range of 
values for which a reduction in the phase variance due to 
squeezing should be observed. 

The last issue which we address in this section is that 
of using arg(C() rather than Bayesian phase estimates. 
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TABLE I: Asymptotic values of cP' sj M j k, under each of the measurement schemes. The numerically estimated values are given 
first, and the analytic predictions are given in brackets. 





adaptive 


heterodyne 


arbitrary squeezing 
squeezing limited to e^^ < 2 
coherent states 


0(0) 
0.351 (l/VS) 
0.497 (1/2) 


0.501 (1/2) 
0.612 

0.705 (l/\/2) 



To estimate the accuracy of the arg(Ct) phase estimates, 
calculations were also performed with Bayesian phase es- 
timates. The system was integrated over a time period 
of lO'^/x, and the phase variance was estimated by av- 
eraging from time 30/x. The functions P{9), x{9) and 
G{9) were estimated by calculating them at 2000 values 
of 9. There is a complication when we take account of 
the variation in the system phase. If the Bayesian anal- 
ysis is performed exactly, the distribution is no longer 
Gaussian in x for given 9. It is not feasible to perform 
the calculation for the full distribution in x and 9, so 
the distribution was approximated by a Gaussian in x. 
The effect of the varying system phase on the distribu- 
tion was approximated by simply adding a spread to the 
phase distribution P{9). 

For each set of data, the phase was estimated both via 
the Bayesian method and as arg(Ct). This allows accu- 
rate comparison of the relative variance. The ratio of 
the variance for the arg(Ct) phase estimates to that for 
the Bayesian phase estimates is shown in Fig. [7| There 
is only a small difference between the two variances; on 
average about 4%, and no more than 9%. In fact, the 
phase estimates obtained via arg(Ct) and the Bayesian 
method are quite close. The mean-square difference be- 
tween these phase estimates is only about 3% of the total 
phase variance. This means that the qualitative results 
obtained for the arg(Cf) phase estimates should also hold 
when Bayesian phase estimates are used. The only dif- 
ference is a reduction in the variance of a few percent. 

The difference between the variances is likely due to the 
fact that the Bayesian estimate uses phase information 
from the noise, whereas the arg(Ct) estimate does not. 
Using the result (|3.53|l we may predict the ratio of the 
variances. This estimate is also shown in Fig. There is 
not exact agreement with the numerically obtained ratio, 
but the estimate is close, particularly for the larger values 

VII. CONCLUSIONS 

We have undertaken a thorough analysis of adaptive 
estimation of a continuously varying phase for narrow- 
band squeezed beams. The problem is characterized by 
photon flux A/", phase diffusion rate k, squeezing band- 
width 7, and degree of squeezing e^^ (i.e. depth of squeez- 
ing of 1 — e^^''). If it is possible to achieve arbitrarily high 
squeezing, then adaptive phase measurements should 
give a variance scaling as (Af / k)~^^^ , an improvement 



> 
o 

■2 0.4 

CO 

0.2 



10° 10= 10^° 

N/k 

FIG. 7: The ratio of the phase variance obtained via arg(Ct) 
to the phase variance obtained via the Bayesian method (solid 
line). The estimated ratio, based on the difference being due 
to the phase information from the noise and using the maxi- 
mum value of (?(r, $ — O), is shown as the dashed line. The 
mean-square difference between the two phase estimates, as 
a ratio to the phase variance for Bayesian estimates, is shown 
as the dotted line. 

over the {M/k)^^^^ scaling for coherent beams. This 
variance is higher than that suggested by the broadband 
analysis in Ref. 0, which had a scaling of {Af / k)^^^^ . 
The reason for this difference is that a time scale of /j 
is required before squeezing is observed. This means that 
the phase information must be obtained from a longer 
time interval than in the broadband case, and the sys- 
tem phase varies by a larger amount over this time in- 
terval. In the broadband case the limiting factor was the 
accuracy of the feedback. 

It is somewhat surprising that the time scale required 
to observe squeezing is /j, rather than I/7, which is 
what we would expect since 7 is the squeezing band- 
width. The rate 7/e'' is the bandwidth of the anti- 
squeezed quadrature, but it is only through the full 
Bayesian analysis that one finds that this is the important 
rate. The difference of is crucial, because the scaling 
of the phase variance would be (A/'/k)~^/^ (as suggested 
by the broadband analysis) without it. The numerical 
calculations verify that the scaling is {Af/K)~^^^ rather 
than (A/'/k)^^/'^, demonstrating that the factor of is 
correct. 

We performed a Bayesian analysis of the phase, for 
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both adaptive and nonadaptive (heterodyne) measure- 
ment schemes. In each case the analysis yields a variance- 
reducing term proportional to E'^, which may be inter- 
preted as the phase information from the coherent am- 
plitude, as well as a term proportional to 7, which cor- 
responds to phase information from the squeezed noise. 
Except in cases where the dominant contribution to the 
photon flux is from the squeezing, the phase information 
from the noise is predicted to be negligible. This means 
that it is possible to obtain accurate phase estimates us- 
ing a simplified method rather than the full Bayesian 
estimate. 

In the experimentally realistic case where there is lim- 
ited squeezing, the reduction in the phase variance is 
approximately e''. Experimentally it is not possible to 
produce beams with arbitrary values of 7, though a lim- 
ited amount of control is possible For larger values 
of A/'/k, the variance is insensitive to the exact value of 
7, and similar variances are obtained for a wide range of 
values of 7. Therefore this issue is not expected to be 
a problem experimentally. A more complete analysis of 
the experimental feasibility of adaptive phase estimation 
on continuous squeezed beams will be given in a future 
work. 



APPENDIX A: EXPECTED VARIANCE 

The standard measure for phase uncertainty, used in 
most of the Refs. S^i^ is the Holevo variance ISj. It 
is defined as |(e*^^)|~^ — 1, where A9 is the error in the 
phase estimate; that is, the estimated phase minus the ac- 
tual system phase. An alternative definition which avoids 
phase estimates with systematic error is Re(e*'^^)^^ — 1 
ml. The Holevo variance under this definition can be 
determined from |lll | 



jAe\ 



P{I\9)e''^d9 



(Al) 



Here / is used to indicate the entire measurement record. 
This may be alternatively written as 



= J p{i)di j p{e\i) 



P{I\Q)dI 



P{9\I)e"^d9 



(A2) 



The measurement record / can be expressed as a func- 
tion of 6 as well as a record of random fluctuations v 
which are independent of @. Therefore we may give the 
expression for (e*'^^) as 



iA0\ 



1-JdeJ Pi.)d. 



P{9\iy,e)e'^d9 



(A3) 



Provided the measurement gives small variance, the only 
values of 9 and Q for which P{i')P{9\i/,Q) is not neg- 
ligible are those for which 9 is close to 9. Thus it is 



reasonable to expand P{9\i', Q) as a series about Q: 

m 

where a, b and c are functions of v and Q. Simplifying 
gives 



iA9\ 



^ J de J F(zy)e-i/(4=)di/ 



(A5) 



Provided the variance of c is small, it is possible to use 
the approximation 

i^'^')^^-^^ J de{{c)er'- (A6) 

If {c)e is independent of 6, (e*'^'') « 1 - l/(4(c)), so the 
variance is approximately l/2(c). 

APPENDIX B: HETERODYNE DERIVATION 

The discretized version of the equations for heterodyne 
detection, H4.5|l and (|4.t)|l is 

Ax + iAy = -7At [a;(l + e)/2 + iy{l - e) /2] + y/2-/Ati)i 

(Bl) 



/ = e'^[y/j/2{x + iy + iE') - (i>i + i>2)/V At]. 
Here the vi, V2 are complex Gaussian random variables 



with {vj) = 



iviVj) but {v*Vj) 



At each time step we update the probability -P(x, y, 9), 
ignoring the constant factor 1/P{I), by 

Pix, y, 9) ^ P{x, y, 9) x P(I\x, y, 9). (B2) 

The probability P{I\x^y,9) is given by 



P{I\x,y,9) cx exp ( 



^^i]2{x + iy + iE') 



exp 



^i^iix - i^r + {y - 



where Ir 



iL, 



'-y) 

(B3) 

3"'*' - iE'. At aU times the 



probability P{x, y, 9) is given by 

P[x,y,9) = P{9)P{x\9)P{y\9), 



(B4) 



where P{x\9) and P{y\9) are normal distributions for x 
and y. P{9) gives the correct probability for 9 averaging 
over x and y. 

We therefore have 



1 



1 



-fAt 



1 



7At 



x^{x + a^47At/2)/(l + al-iAtl2), 

y^{y + ally^At/2)/{l + al^Atl2). (B5) 
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In the infinitesimal limit, 

dcrl = -cjtldt/2, dal = -a^'ydt/2, 
dx = al{I^ - x)'ydt/2, dy = a^Iy - y)'ydt/2. (B6) 

For the change in the probability distribution for 9 we 
have, in the infinitesimal fimit 

P(2) (6) a P(0) exp I - ^ [(4 - + (/, - I . 

(B7) 

Now wc take account of the increments in x, y and 6. 
We have the restriction that 



Therefore the steady state values are 



j>i + j>2 = VAt[V7/2(x + iy + iE') - Ie-'% (B8) 

The mean of Di is therefore half of the RHS, and the 
second-moment is 1/2, rather than 1. Therefore, the dif- 
ference equation may be expressed as 



Ax+iAy = -[{x-iy)e+I^ + iIy]jAt/2 + i>./^, (B9) 

where /> is a complex Gaussian random variable. Consid- 
ering the deterministic part of the increment first, and 
taking the infinitesimal limit 

dal = -(^Indt, dal = (^P^dt, 

dx = —{Ix + ex)^dt/2, dy = —{ly — ey)^dt/2. 

(BIO) 

To take account of the stochastic increments, we simply 
add the appropriate variances: 

dal = ldtl2, dal = ^dtl2. (Bll) 

Overall we have the increments 

dal = [1 - 2£a2 - a^]7rfi/2, 

dal = [1 + '^^^l - ^t]^dt/2. (B12) 



ct2 = v'e^ + 1 - £, (7^ = Ve^+T -h £. (B13) 



The total increments in x and y are 



dx = -[x{al +e) + {l- al)Ix]^dt/2, 
dy = -Wl -£) + (!- al)Iy]^dt/2. 



(B14) 



The solutions are 



4e 



dv, 



x = {al-l){j/2) I 
Jo 

rt 

y=(a2-l)(7/2) / 7ye(T/2)v^(''-«)d^;. (B15) 
Jo 



Now (/) = ie^sfij2E\ so {h+Hy) = i(e'(®-^) - l)E' , 
and 



{x) = sm{6 - e)E 



.crl-1 



(y) = [cos{0 - e) - 1]^;' 

The solutions for x and y are 



(B16) 



x{t) = ^27 /* 6-^(1+") ("-*)/2Re(j/i)du, 
Jo 

y{t) = ^ [ eT(i-=)("-*)/2lni(!/i)du. (B17) 
Jo 



In terms of these quantities. 



= x{t) cos{e -e) + {y{t) + E') sm{6 - 6) - ^/2Pi{Rc{vi + V2) cos((9 Im(!/i -|- V2) sin(6l - 6)} 

Iy{t) = -x{t) sm{e -Q) + {y{t) + E') cos{e - Q) - E' - ^/2ff{lm{ul + V2) cos{e - 6) - Re{ui + U2) sin(6l - 9)}. 

(B18) 

Therefore 



{Ix{u)Ix{v)) 

{Iy{u)Iy{v)) 



{{x{u)x (?;)) - ^/2/^{x{u)Rc[iyi{v)] + x{v)Rc[iyi{u)])} cos^{0 - e) 

+ {{yiu)yiv)) - y^{yiu)lm[iyiiv)] + y(w)Im[j/i(u)]) + {E' f} sin^O - 9) + (2/7)(5(« - v) 
{{x{u)x{v)) - y/2P'f{x{u)Rc[i^i{v)] + x{v)Re[iyi{u)])} sm^{e - 9) 
+ {{yiu)yiv)) - ^^^2p}{yiu)lm[iyi{v)] + y(w)Im[:/i('«)])} cos^(6' - 9) 



+ {E'y[cos{d - 9) - 1]^ + (2/7)(5(u - v). 



(B19) 
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Substituting gives 



{Ix{u)Ix{v)) = e ''e 



- e) + {e'-e-''(i-^)l"-"l/2 + [E'f] sin2((? - 6) + {2h)5{u - v) 
+ (£;')^[cos(6i - 9) - 1]^ + {2/^)5{u - v). (B20) 



/\2r //I ii2 



{Iy{u)Iy{v)) ^ e ''e 'v-'-^i- "I'-sin' 



Using this result gives 
{{h - x)^ + (ly - y?) = const. + {E'f ( 1 



2 i\2 /'„2 i\2 



9) + 1 



a^-l 



a2 + 1 + 1 



[008(6* - 9) - 1]2 

(B21) 



Here the terms that do not depend on or 9 have been coUected into the constant. Expanding to second order in 



— 9) and using the equihbrium values of cr^ and cr^ gives 



^((4 ~ xf + (/, - yf) = const. + 



2r ' 



ef + 2^dt coshr- (e^ + i)-i/2 
I 



-Qf 



(B22) 
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